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ABSTRACT 

Sand bed roughness was measured in the nearshore during a variety of 
hydrodynamic conditions using an array of seven 1 MHZ sonar altimeters mounted on the 
WESP (an amphibious vehicle used for measuring large-scale bathymetry) during the 
COAST 3D Experiment 1998, at Egmond aan Zee, The Netherlands. Corollary waves 
were modeled using the Thronton and Guza (1983) wave transformation model. Wave 
height and current measurements were made in the surf zone using pressure sensors and 
electromagnetic current meters, and offshore wave heights were measured using an 
advectional wave-rider buoy. Measurements of sand bed roughness showed patterns 
similar to those observed by Clifton et al. (1971), but are highly variable both spatially 
and temporally with dependence on large scale morphology and wave and current 
conditions. Mobility number ('F) calculated from the modeled wave field and also from 
measured currents revealed that roughness is a function of Y. Roughness was observed 
to be highly variable at low 'F (< 100) calculated from the modeled wave heights. As 
increased (100-150), roughness decreased gradually, but was still highly variable. As T 
reached values greater than 150, roughness was restricted to the lowest observed values 
(< 3 cm) implying planar beds under sheet flow conditions. 
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I. INTRODUCTON 



Barred beaches are composed of highly variable and complex bedforms (Clifton 
et al., 1971). The formation and distribution (spatial and temporal) of these bedforms is 
poorly understood, and a better understanding would be a valuable tool for determining 
bottom friction and resulting wave energy dissipation, current generation in the nearshore, 
sediment transport and resulting bathymetry changes, as well as for interpretating paleo- 
environments. 

Significant information is available on small-scale morphology in steady 
unidirectional currents, e.g. rivers and estuaries (Middleton and Southward, 1984, 
Fredsoe and Deigaard, 1992). In general, bed state is separated into different regimes 
(e.g., ripples, flatbed, etc.) which are a function of flow strength, flow depth and grain 
size. Similarly, Clifton et al. (1971) described the spacially varying bedforms on a non- 
barred beach in terms of flow regimes (Fig. 1). They observed that well outside the surf 
zone the bed is dominated by straight asymmetric ripples, 2-5 cm high, with 10-20 cm 
wavelengths. Immediately seaward of the breaker line in the zone of wave build-up, they 
described the occurrence of landward-oriented lunate megaripples (30-100 cm high, 1-5 
m long). Beneath the breakers they observed a planar bed. The inner portion of the surf 
zone is characterized by large-scale, long-crested, symmetric ripples described as having 
heights of 15-20 cm and lengths of 30-60 cm, as well as large megaripples. Planar beds 
characterized the swash zone. Clifton et al. (1971) suggested that this landward sequence 
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is the result of an increase in flow regime associated with a shoreward increase in orbital 
velocity and velocity asymmetry at the bottom. 

Clifton (1976) developed a conceptual bedform model based on maximum bottom 
orbital velocity, velocity asymmetry, median grain size and wave period, which consists 
of four main regimes including no sediment flow, symmetric bedforms, asymmetric 
bedforms and flat bed (sheet flow). Davidson-Amott and Greenwood (1974) also studied 
bedforms on a barred beach and developed a basic regime model which divides each bar 
and trough into four subregions, and corresponding flow regimes analogous to Clifton 
(1976). 

These studies employed SCUBA to observe bedforms in situ and therefore were 
conducted during relative calm conditions. During storms, the turbulent environment of 
the nearshore is impassable, thus observations during high energy conditions were 
impossible. More recently, acoustics have been used. For example, stationary sonar 
altimeters have been used to investigate bedforms in the nearshore (Dingier and Clifton, 
1984, Hay and Wilson, 1993) and to measure accretion and erosion in the surf zone 
(Gallagher et al., 1998) and in intermediate water depth (8 m) (Wright et al., 1986). 
Dingier and Inman (1976), Dingier and Clifton (1984), and Greenwood et al. (1993) used 
a track-mounted vertical sounder to obtain 2-D profiles of bed elevation at a single 
location inside the surf zone. Hay and Wilson (1994) used a rotary sidescan sonar to 
document the continuous evolution of bedforms through time over a 10 m^ area on the 
crest of a nearshore bar where they observed the occurrence of ripples, cross ripples, 
meganpples, and transition between ripple types during the waning stages of the storm. 



2 



Thornton et al. (1997) measured small-scale morphology related to wave and 
current parameters across the surf zone using a single altimeter mounted on an 
amphibious surveying vehicle. Gallagher et al. (2000) made similar measurements with 
an array of seven sonar altimeters mounted on the same moving platform to characterize 
the sand bed roughness in the nearshore over a 500 m X 700 m area. Here similar 
measurements of bedforms for a different beach are presented. In Chapter II the 
experiment is described, and in Chapter El analysis techniques will be documented. In 
Chapter IV, the observations will be presented, and in Chapter V they will be discussed. 
The objective of this study is to characterize sand bed roughness and its spatial and 
temporal variability in the nearshore, and to quantify the degree of roughness by 
comparing to wave conditions. 
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II. COAST 3D EXPERIMENT 



The data presented were obtained as part of COAST 3D, a comprehensive 
nearshore field experiment (Mulder, 1997) conducted over a six- week period in October- 
November 1998 at Egmond aan Zee, the Netherlands (Fig. 2). The objective of the study 
described here was to measure the spatial and temporal distribution and variability of 
bedforms in the nearshore under a wide variety of conditions. 

The WESP, a 15 m high, amphibious, motorized, three- wheeled vehicle made 
daily surveys of the large-scale bathymetry (e.g. sand bars) using DGPS (Fig. 3). An 
example of a two-dimensional bathemetric contour map produced from the data collected 
by the WESP describing the large-scale morphology is shown in Figure 4. The vertical 
and horizontal error of the WESP bathemetric surveys is estimated at 10 cm (Ruessink et 
al., submitted), and the alongshore separation of the survey lines was about 50 m. 
Bathymetry measured by the WESP was sampled at approximately 1 Hz, (although the 
data were not evenly spaced in time). The spatial separation of data is dependent on the 
speed of the WESP. The WESP speed is 0(1 m/s) resulting in a cross-shore resolution of 
0(1 m). 

The elevation of the small-scale morphology, which the WESP cannot resolve, 
was measured using an array of 7, 1 MHz sonar altimeters attached to the rear of the 
WESP (Fig. 5) following the methods of Gallagher et al (1996), Thornton et al. (1997), 
and Gallagher et al. (2000). The altimeters were sampled at 48 Hz so the cross-shore 
resolution (also dependent on WESP speed) is 0(2 cm). The altimeters have a much 
higher vertical resolution 0(mm) than does the WESP DGPS, and the closest spacing 
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between altimeters is 15 cm. The data were corrected for pitch, roll and yaw of the 
WESP using measured heading and tilts. The cross-shore profiles produced using the 
altimeters are closely spaced with vertical resolution of about 3 cm (Gallagher et al., 
1998) and cross-shore resolution of about 5 cm and alongshore resolution of about 40 cm. 
Thus altimeters are capable of measuring the small-scale bedforms (i.e., large ripples and 
megaripples). 

The field site forms part of the 120 km long uninterrupted Holland coast and faces 
the North Sea. The nearshore subtidal morphology is characterized by a shore-parallel 
double sandbar system with a mean slope of 1:100 (Ruessink et al., 2000). Sediments at 
the field experiment near Egmond are well sorted and composed of fine to medium sand 
with a mean grain size between 0.25 and 0.35 mm (Krone et al., 1997). The inner bar 
often contains significant quasi-periodic alongshore variations with a wavelength varying 
between 350 and 900 m (Ruessink et al., 2000). The outer bar can be described as being 
either straight or containing rhythmic features with length scales of more than 1 km 
(Ruessink et al., 2000). 

Offshore wave conditions (Tablel, Fig. 6) were acquired from an advectional 
wave rider buoy located approximately 5 km offshore in 16 m water depth. The first 14 
days of the COAST 3D field experiment were characterized by an increase in significant 
wave height, H1/3, from 1 to 3 m (Fig. 6a). During the next six days, large waves 
(H]/3 = 2-5.5 m) were observed. Afterwards, the offshore wave conditions decreased 
considerably with the exception of day 26. Significant wave period, T1/3, ranged from 3 
to 10 s (Fig 6b), mimicking the development of the H1/3. Wave directions during the 
experiment were from the southwest (-45®) or from the northwest (445®) (Fig. 6c). The 
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tide is semi-diurnal with a mean range of 1.5 m. Storm surge levels periodically raised 
water levels by more than 1 m causing low-water levels to be higher than ordinary flood 
levels (Ruessink et al., 2000). Figure 6d is offshore water level measure in NAP, Dutch 
ordnance level. 



7 



TfflS PAGE rNTENTIONALLY LEFT BLANK 



8 



III. METHODOLOGY 



The WESP survey and altimeter measurements were combined to obtain highly 
resolved bottom profiles at the Egmond field site during the COAST 3D Experiment 
(Fig. 7). GPS data were linearly interpolated to 48 Hz and a median filter was used to 
remove bad returns from the altimeter data (Thornton et al., 1997). After correcting 
altimeters for pitch, roll and yaw, profiles of cross-shore distance versus WESP survey 
depth plus height from the bed to the altimeters were obtained (Fig 7a). Cross-shore 
profiles of bedforms only were generated by subtracting smoothed WESP profiles of the 
large scale bathymetry. 

This study will focus on the root mean square (RMS) roughness of the seafloor 
after the method of Gallagher et al. (2000), although the array is capable of resolving 
cross- and alongshore wavelengths of bedforms. Cross-shore profiles RMS of roughness 
are calculated using demeaned, 80 % overlapping, 25 m-long sections of the bedform 
profile. For each survey grid line, the seven roughness profiles (one from each sonar) are 
averaged to give a single roughness profile (Fig. 7). Thus, for each daily WESP survey, 
which consisted of 20 grid lines (10 grid lines for a small survey), a map of RMS 
roughness is generated (see for example Fig. 8). In Figure 8, the map of roughness is 
overlain by contours of large-scale bathymetry measured by the WESP. 

Bad RMS roughness values were removed from the profiles prior to mapping. 
Roughness points were determined to be bad if a particular point was not consistent with 
surrounding points, e.g., if a point was > 1.5 time the surrounding points it was 
considered anomalous (this threshold was found by trial and error). Also, roughness 
points were sometimes anomalous owing to aliasing of rapid changes is large-scale 
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morphology. For example, a steep slope on a sand bar was poorly resolved by the 
WESP, and not properly filtered as a large-scale feature. Thus the bar was included in 
estimates of bedform roughness. (Better filtering techniques could correct this aliasing 
problem.) Each roughness profile was compared to its corresponding bedform-only 
profile and cross-shore depth profile, and any bad points were removed. 
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IV. OBSERVATION 



The RMS roughness (here on referred to as "roughness") patterns in the nearshore 
are overlain by bathymetry contours in Figures 8-14. The typical pattern (e.g. Fig. 11, 
Hi/3=0.98 m) consists of roughness on the shoreward slope of the outer bar with the outer 
trough characterized by less roughness. Roughness increases again on the seaward face 
of the inner bar and is greater than that observed along the shoreward face of the outer 
bar. Roughness decreases near the crest of the inner bar and increases again shoreward of 
the inner bar. 

However, observed patterns of roughness in the shallow water are highly variable. 
For example, a short survey (500 m by 300 m. Fig. 12, Hi/3=0.7 m) revealed significant 
roughness on the crest of the inner bar and a relatively smooth bed on the seaward slope 
of the bar. Although this survey does not cover the whole area, the variation from the 
above pattern can be seen. In addition, breaks in the inner bar, possibly owing to rip 
currents, show areas of anomalously high roughness (Fig. 11, Fig. 13, and Fig. 14). High 
roughness in rip current channels has been observed by Thornton et al (1997) and Vincent 
and Osbourne (1993). 

The typical roughness patterns observed during this study are analogous to 
Clifton’s (1976) observations (Fig. 1). His conceptual model shows the transitions of 
small-scale morphology from wave ripples to megaripples on the seaward slope of the 
bar, to a planar bar crest due to the increase in wave intensity, and to wave or megaripples 
again in the trough where wave intensity decreases after waves have crossed the bars. 
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This general sequence of bedforms with increasing flow strength described by Clifton 
(1976), and observed by others (e.g. Davidson-Amott and Greenwood, 1976, Hay and 
Wilson, 1994, Thornton et al., 1997, and Gallagher, 2000) is: 

no movement - ripples (small) - megaripples - plane bed. 

This sequence is typical for fine sand (<0.2 mm) (Davidson-Amott and Greenwood, 
1976). Davidson-Amott and Greenwood (1976) observed that as grain diameter increases, 
the occurance of megaripples increases, and that of ripples decreases until coarse sand 
(1.0 mm and greater) when ripples are no longer formed and the sequence becomes: 
no movement - plane bed - megaripples - plane bed. 

There have been variations reported on these typical bedform sequences. For example. 
Southard and Boguchwal (1973) have shown that there is a transitional area between fine 
and coarse sand in which the succession of bedforms is more complex: 

no movement - ripples - plane bed - megaripples - upper flat bed. 

The bedform sequence observed at Egmond, where median grain size is 0.3 mm, 
appears to be typical for fine sand. For example, typical cross-shore transects are found 
in Figures 15-19. (Example transects are taken from alongshore location y=775 m 
because they are co-located with the pressure and current measurements, see section V.) 
On these days, roughness is low (<5 cm) on the crest of the outer bar where waves are 
shoaling, and possibly breaking, suggesting planar bed under sheet flow conditions (Fig. 
15d,e and 19d,e). As water depth increases over the shoreward face of the outer bar and 
trough, roughness increases (5-12 cm) but decreases again (< 5 cm) in the outer trough. 
An increase in roughness (7-10 cm) is seen as water depth decreases again on the seaward 
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face of the inner bar. Along the crest of the inner bar, a decrease in roughness is seen 
owing to increased flow strength as waves break, perhaps generating sheet flow or small 
ripples. The inner trough is characterized by large, but highly variable roughness, what 
Clifton (1976) called the “inner rough” regime, probably owing to the complex flow field 
in the inner surf zone (Fig. 16d and 19d). 

While the above examples show striking similarities to the Clifton (1976) 
bedform distribution sequence, there are variations that deserve attention. For example, 
during the highest waves of the experiment, (Table 1, Fig. 18c) the crests of the bars are 
planed off as well as the bar slopes, and in the troughs the roughness is very high 
(8-12 cm) indicating megaripples. There appears to be little or no transition down to 
what may be analogous to small ripples associated with a low flow regime. 

Another variant on the bedform sequence conceptual model (Clifton, 1976) is the 
effect of steady currents (longshore, tidal, rips) on roughness. Between 24 October and 8 
November, there were a series of storms, which cut rip current channels through the inner 
bar (Fig. 13 and Fig. 14). Roughness in the rip channels (Fig. 13, y=775 m and Fig. 14, 
y=700 m) is noticeably higher (5-8 cm) compared to the roughness (< 5 cm) on the 
adjacent bar crests to the north and south. The dominant process controlling bedform 
generation and movement in the rips is most likely the unidirectional current, and wave 
oscillatory currents are of lesser significance (Davidson- Amott and Greenwood, 1976). 
Bedforms generated by combined wave-current flows occur in many coastal sub- 
environments, and Nielson (1992) observed that in shallow bar troughs, or rip feeder 
channels, conditions are cuirent-dominated, with bedforms resembling current dunes. 
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Thornton et al. (1997) observed seaward facing lunate and straight-crested megaripples at 
Duck, NC suggesting that the offshore rip current dominated bedform processes over the 
oscillatory wave-generated flow. 



14 



V. Discussion 



Clifton (1976) hypothesized that different bed form regimes are caused by 
variations in the velocity field as a function of wave transformation and breaking across 
the surf zone. To investigate the relationship between the roughness and the wave field, 
wave height changes were modeled across all survey lines using the model of Thornton 
and Guza (1983). The wave height transformation model was initialized using input of 
hourly averaged RMS wave height (Hrms), peak frequency and mean direction measured 
at a wave rider buoy in 16 m of water. The model has two free parameters, B and Yb 
corresponding to the intensity of wave breaking and breaking waves criterion, which were 
calibrated using Hnns values measured by the array of pressure sensors across the inner bar 
(Fig. 15a- 19a) and agree well with values found in the literature (Thornton and Guza, 
1983). 

Predicted Hrms versus cross-shore distance (at alongshore location y=775 m) are 
shown in Figures 15a- 19a, and the comparison with measured Hrms (asterisks) is good. 
The predicted cross-shore values of Y=Hrms/h (h=depth) are shown in Figures 15b- 19b to 
vary with the large-scale bathymetry (Fig. 15e-19e). It is hypothesized that the bed 
roughness is related to this normalized wave height. RMS roughness is plotted versus y 
for all values of the COAST 3D experiment in Figure 20, and a clear relationship between 
Y and roughness is not evident. 
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It was further hypothesized that mobility number (T), a dimensionless ratio of 
fluid forces (lift and drag) to gravitational force on sediment particles under waves, is a 
function of roughness. Nielsen (1992) summarized 'F for waves alone, 

(s - l)gd 

where Urmsb is the rms wave induced velocity at the bottom, s is the specific gravity of 
quartz, g is acceleration due to gravity, and d is grain diameter (0.3 mm). The RMS wave 
induced bottom velocity is calculated from linear theory 

* 2 sinh kh 

where Hrms is the wave height calculated from the Thornton and Guza (1983) model, and 
k is wavenumber. It is pointed out that the V2 is due to relating the Hrms— '^Surms as 
calculated assuming a Rayleigh wave height distribution. Wavenumber is related to 
radial wave frequency (to) using the dispersion relationship from linear wave theory, 

£y‘ = gk tanh kh. 

Mobility number calculated in this way is shown in Figures 15c- 19c and, like y, shows a 
strong relationship with the large-scale morphology. The relationship between T and 
RMS roughness is shown in Figure 21. As found by Dingier and Inman (1976), a 
decrease ripple height (roughness) corresponds to higher T. However at lower T the bed 
roughness is highly variable. 
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Sediment movement is initialized when the critical shields parameter (dc) is 



exceeded. The Shield’s parameter (6) is given by 






where is an empirical wave friction factor (/iv=0.1 is a conservative estimate), and 'F is 
mobility number as defined above. The critical Shield’s parameter {Be) was calculated to 
be 0.4 in this study using the method of Komar and Miller (1973). Based on our 
calculations, Be was always exceeded during this experiment. 

Nielsen (1981) summarized the results of wave ripple measurements in the field 
(Inman, 1957; Dingier, 1974; Miller and Komar, 1980) relating ripple wavelength (X) and 
wave orbital excursion amplitude, a, to mobility number. Nielsen (1981) suggested for 
values of 'F greater than 200, \la becomes minimal, implying that ripples are planed off. 
Dingier (1974) suggested, based on field observations, that ripples are planed off when 'F 
reaches a value of 240, when considering sediments of almost entirely quartz sand. 
Nielsen (1981) also found that where sediments have densities close to that of quartz, it is 
reasonable to expect relative ripple height (r)/a) to be a function of 'F as well. Nielsen 
(1981) suggested that relative ripple height goes to zero for T>156. 

In this study, as T increased to its highest values (>150), roughness is limited to 
the lowest observed values (< 3 cm. Fig. 21), which generally agrees with the thresholds 
for planar bed proposed by Dingier (1974) and Nielsen (1981). The roughness is highly 
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variable throughout a wide range of 'P, and in general the roughness decreases with 
increasing 'F. The highest roughness values are typically observed where 'P<100. 

Flume studies of unidirectional flow have shown that the transition from small 
ripples to large ripples (megaripples) is complete over a narrow range of only a few 
centimeters per second in flow velocity. Then, with increasing velocity, megaripples are 
observed to become lower and more rounded over a wide interval of flow velocities 
(Middleton and Southard, 1984). While it is difficult to compare unidirectional flow 
studies (such as rivers and estuaries) to the complex flow regime of the nearshore, 
perhaps the tapering off of roughness for ^>100 is analogous to the transition from 
megaripples to planar beds. However, no sharp transition is observed between ripples 
and megaripples. 

Predicted 'F is based on cross-shore oscillatory velocity only, thus any 
contribution to sediment mobility owing to steady currents is not considered. Using 
measured velocities from electromagnetic current meters (co-located with pressure 
sensors at asterisks on Fig. 15e-19e), total 'F is calculated 

— 9 

.j, _ 

(s-l)gd 

where u = (urms + U) + V with U and V are the mean cross-shore and alongshore 
currents and Urms is the square root of the cross-shore velocity variance. Alongshore 
velocity variance is assumed negligible. 

The current meter-derived T values are plotted against roughness measured in the 
vicinity of the sensors (Fig. 22). While there are only few points, a noticeable trend 
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exists, and Clifton’s (1976) conceptual model is supported: small roughness elements 
(wave ripples) exist for small values of T, larger bedforms exist for intermediate 'F 
values. According to Nielsen (1992) and Dingier (1974), a reduction of roughness would 
be expected above about 'F=200 corresponding to the transition to sheet flow. 
Unfortunately, no data exist for these conditions. 

A mean roughness value was obtained for each day by averaging all roughness 
values for each day in both the along- and cross-shore direction. The same cross-shore 
extent was averaged for each survey, but alongshore extent was allowed to vary assuming 
alongshore homogeneity and cross-shore inhomogeneity. Mean roughness is plotted 
versus significant wave height for each day in Figure 23. Similar to purely steady flow 
bedform regime studies, for low energy and high energy conditions overall roughness is 
lower, corresponding to rippled beds or sheet flow, respectively. However for lowest 
wave conditions roughness variability is the highest, which is consistent with data in 
Figure 21. Although confidence in these average roughness values in high (Table 1), 
possible uncertainty (3 cm, see Chapter II) associated with these measurements makes 
this observation barely significant. 

Grain size and grain size distribution, although important in sediment mobility, 
were not considered here. In most coastal environments, the general grain size 
distribution from offshore to onshore is fine to coarse, however grain size distribution can 
also change with a change in wave conditions. Because grain size affects the size or 
elimination of bedform regimes (Davidon-Amott and Greenwood, 1976), including 
spatial and temporal distribution of sediment may be important. Grain size variation 
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needs to be considered in future studies of bedforms in the nearshore to help explain the 
observed complexity. 

In addition, the assumption that the bedforms are in equilibrium with the present 
conditions may not be valid. Although relict ripples are unlikely in the high energy 
environment of the nearshore, there may be a lag between changing flow conditions and 
bedform development. Thus, a lag between forcing and morphology changes should also 
be considered in future studies. 
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VI. CONCLUSIONS 



Sand bed roughness in the nearshore was measured in a variety of hydrodynamic 
conditions on a double-barred beach during the COAST 3D experiment at Egmond aan 
Zee, The Netherlands. The root-mean-squared (RMS) roughness was measured using an 
array of seven sonar altimeters mounted on the WESP. Corollary wave conditions were 
modeled using the Thornton and Guza (1983) wave transformation model, and pressure 
sensors and electromagnetic current meters in the surf zone acquired wave height and 
current measurements. 

Sand bed roughness patterns observed in the nearshore at Egmond exhibited 
patterns similar to those described by Clifton et al. (1971), but were found to be highly 
variable (spatially and temporally). Roughness appears to be dependent on the large-scale 
morphology and influenced by the hydrodynamic conditions. 

Mobility number was calculated from bed velocities derived from modeled 
wave heights using linear theory. For T > 150, roughness was restricted to the lowest 
values (< 3 cm) observed during the experiment, implying the occurrence of plane beds 
under sheet flow conditions. These values of T for when bed forms are planed off agree 
with those found in the literature. However, for conditions of low 'P, roughness was 
shown to be highly variable. This variability could be owing to strong steady currents 
(rips, alongshore, tidal, etc.), which were not accounted for in the modeled mobility 
numbers. 

Current velocity was measured with four current meters located on the inner bar. 
Mobility numbers calculated from the measured current velocity show that roughness has 
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a strong dependence on the current-derived 'F. During periods of low 'F, roughness was 
observed to be 3-5 cm, probably owing to vortex ripples. As 'F increased, roughness 
increased to 5-15 cm owing to large ripples and megaripples. No data were available for 
'F > 200, thus the expected decrease in roughness was not observed. 

Other factors are known to influence the development of sand bed roughness, but 
were not considered in this study. Variations in sediment composition and distribution is 
thought to influence sand bed features in terms of ripple wavelength, slope, and extent of 
ripple types. Interaction of currents, such as longshore, rip, and tidal have been shown to 
dominate bed form processes over that of oscillatory currents in certain areas of the 
nearshore. Remnant roughness, or a time lag between changing flow conditions and bed 
form development, were not considered during this study. Further investigation into 
these factors will contribute to the understanding of sand bed processes in the nearshore. 
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Day 


Significant Wave 
Height (m) 


Significant Wave 
Period 
(s) 


Mean 

Roughness 

(cm) 


Roughness 
Standard 
Deviation (cm) 


95 % Confidence 
Interval (cm) 
(±) 


19 Oct 


1.8 -2.0 


6.3 - 6.4 


6.0 


1.7 


0.05 


22 Oct 


2.3 - 2.4 


7.2 


5.7 


2.7 


0.08 


24 Oct 


4.7 - 5.2 


8.6 - 9.0 


4.8 


0.8 


0.01 


06 Nov 


0.7 


4.8 - 5.2 


- 


- 


- 


08 Nov 


0.9- 1.4 


3.8 - 5.7 


5.2 


0.9 


0.03 


15 Nov 


0.7- 1.0 


5.8 -6.5 


3.5 


0.6 


0.01 



Table 1. Offshore wave conditions and mean roughness for COAST 3D roughness 
surveys. 
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APPENDIX B. FIGURES 
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Figure 1. Distribution of small-scale morphology observed on a barred beach (from 
Clifton, 1976). 
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Figure 2. Map of the COAST 3D field site, Egmond aan Zee, The Netherlands (from 
Ruessink, 2000). 
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Figure 3. WESP (Water en strand Profiler) resting on the beach during COAST 3D 
Experiment. The WESP is about 14 m tall and 8 m wide at the rear wheels (forward is to 
the right), and the diameter of the tires is about 1.7 m. 
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Figure 4. Example of a two-dimensional contour map produced from the data collected 
by the WESP on 15 November 1998. Asterisks denote the location of electromagnetic 
current meters and pressure sensors. 
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Figure 5. Close up of the array of seven sonar altimeters mounted on the rear cross 
member of the WESP (Fig. 3). The sonar transducers are about 70 cm above the sand. 
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Figure 6. Offshore wave conditions during the COAST 3D Experiment, a) significant 
wave height, b) significant wave period, c) wave direction, d) offshore water level, 
where NAP is Dutch ordnance level. 
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Figure 7. a) Example of a cross-shore profile from a single sonar altimeter with elevation 
below mean sea level vs. cross-shore distance, b) Example of 25 m-long piece of Fig. 7a. 
These data are detrended, and the root-mean-squared is calculated to give RMS 
roughness, which for this section of bedforms, is 3 cm. The RMS is calculated from 
overlapping sections (as in 7b), to produce this RMS roughness profile from the depth 
profile in Fig. 7a. 
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Figure 8. RMS roughness map overlain by bathymetry from data collected by the WESP 
for 15 Nov at the COAST 3D field site. The beach is on the right-hand side of the plot. 
Asterisks represent the location of current meters and pressure sensors. The cross-shore 
profile in Fig. 7 is take from y=500 m on this map. 
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Figure 9. RMS roughness map overlain by bathymetry from data collected by the WESP 
for 19 Oct. at the COAST 3D field site. The beach is on the right-hand side of the plot. 
Asterisks denote location of current meters and pressure sensors. 
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Figure 10. RMS roughness map overlain by bathymetry from data collected by the 
WESP for 22 Oct. at the COAST 3D field site. The beach is on the right-hand side of the 
plot. Asterisks denote location of current meters and pressure sensors. 
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Figure 11. RMS roughness map overlain by bathymetry from data collected by the 
WESP for 24 Oct. at the COAST 3D field site. The beach is on the right-hand side of the 
plot. Asterisks denote location of current meters and pressure sensors. 
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Figure 12. RMS roughness map overlain by bathymetry from data collected by the 
WESP for 6 Nov. at the COAST 3D field site. The beach is on the right-hand side of the 
plot. Asterisks denote location of current meters and pressure sensors. 
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Figure 13. RMS roughness map overlain by bathymetry from data collected by the 
WESP for 8 Nov. at the COAST 3D field site. The beach is on the right-hand side of the 
plot. Asterisks denote location of current meters and pressure sensors. 
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Figure 14. RMS roughness map overlain by bathymetry from data collected by the 
WESP for 15 Nov. at the COAST 3D field site. The beach is on the right-hand side of 
the plot. Asterisks denote location of current meters and pressure sensors. 
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Figure 15. Cross-shore profiles for 19 Oct (at y=775 m on Fig. 9). a) modeled wave 
height; asterisks denote pressure meter measured wave heights, b) modeled gamma; 
c) mobility number calculated from modeled wave heights, d) roughness profile; and 
e) cross-shore depth profile, where asterisks denote location of pressure sensors and 
current meters. 
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Figure 16. Cross-shore profiles for 22 Oct (at y=775 m on Fig. 10). a) modeled wave 
height, asterisks denote pressure meter measured wave heights, b) modeled gamma; 
c) mobility number calculated from modeled wave heights, d) roughness profile; and 
e) cross-shore depth profile, where asterisks denote location of pressure sensors and 
current meters. 
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Figure 17. Cross-shore profiles for 24 Oct (at y=775 m on Fig. 11). a) modeled wave 
height, asterisks denote pressure meter measured wave heights, b) modeled gamma; 
c) mobility number calculated from modeled wave heights, d) roughness profile; and 
e) cross-shore depth profile, where asterisks denote location of pressure sensors and 
current meters. 
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Figure 18. Cross-shore profiles for 8 Nov (at y=775 m on Fig. 13). a) modeled wave 
height, asterisks denote pressure meter measured wave heights, b) modeled gamma; 
c) mobility number calculated from modeled wave heights, d) roughness profile; and 
e) cross-shore depth profile, where asterisks denote location of pressure sensors and 
current meters. 
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Figure 19. Cross-shore profiles for 15 Nov (at y=775 m on Fig. 14). a) modeled wave 
height, asterisks denote pressure meter measured wave heights, b) modeled gamma; 
c) mobility number calculated from modeled wave heights, d) roughness profile; and 
e) cross-shore depth profile, where asterisks denote location of pressure sensors and 
current meters. 
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Figure 20. Roughness versus gamma for the whole COAST 3D Experiment. 
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Figure 21. RMS roughness versus mobility number derived from predicted wave heights. 
The plane bed thresholds of Dingier (1974), and Nielsen (1981) using relative ripple 
height and relative ripple wave length are shown. The curve illustrates gradual decrease 
in roughness as mobility number increases to plane bed threshold. Dashed line is the 
lower limit at which the sonar altimeters can resolve ripples during this study. 
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Figure 22. RMS Roughness versus mobility number derived from current meter 
measurements. 
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Figure 23. Mean (over a whole survey) RMS roughness versus significant wave height. 
Bars represent standard deviation of RMS roughness. 
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